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(N i ABSTRACT 

A photometric method for quantifying deviations from axisymmetry in optical 
images of disk galaxies is applied to a sample of 32 face-on and nearly face-on spirals. 
The method involves comparing the relative fluxes contained within trapezoidal sectors 
arranged symmetrically about the galaxy center of light, excluding the bulge and/or 
barred regions. Such a method has several advantages over others, especially when 



quantifying asymmetry in flocculent galaxies. Specifically, the averaging of large 
regions improves the signal-to-noise in the measurements; the method is not strongly 
affected by the presence of spiral arms; and it identifies the kinds of asymmetry that 
are likely to be dynamically important. Application of this "method of sectors" to 
i?-band images of 32 disk galaxies indicates that about 30% of spirals show deviations 
from axisymmetry at the 5a level. 

Subject headings: galaxies: photometry — galaxies: structure — techniques: 
photometric 



1. Introduction 

Traditionally, studies of spiral galaxies have concentrated on understanding the dynamics 
of axisymmetric disks with simple spiral structure. The high contrast, well-delineated spiral 
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structure in ideal "grand design" spiral galaxies is fairly well understood as arising from tidal 
effects. However, as Baldwin et al. (1980) (hereafter BLS) first pointed out, ideal (that is, 
symmetric) galaxies outside of the grand design class are more uncommon than generally believed. 
It is now becoming established that perfect symmetry is not nearly as all-pervasive in either 
visible galaxies or in H I disks as previously assumed. 

Asymmetries are generally interpreted (see, for example, Junqueira & Combes 1996) as the 
excitation and superposition of m = 1 spiral modes over the dominant m = 2 mode, or as local 
phase-shifting (with radius) of the m = 2 mode initiated by gravitational perturbations due to 
companion galaxies, even when interacting galaxies are not in evidence. This idea can be extended 
to include those interactions which result in minor mergers, or gas accretion from the surrounding 
intergalactic medium, allowing limits to be placed on the total galaxy accretion rate (Zaritsky & 
Rix 1997, hereafter ZR). 

Alternatively, asymmetry in the H I gas velocity field has been interpreted as an indicator 
of more subtle dynamics occurring within the galaxy itself. It is beginning to be understood 
that stationary solutions of the hydrodynamics of baryonic and dark matter components may be 
non-axisymmetric. Shoenmakers et al. (1997), for instance, interpret optical asymmetry as an 
indicator of asymmetry in the overall galactic potential, and therefore an indicator of the (triaxial) 
spatial distribution of the dark matter in a galaxy. Jog (1997) has studied the orbits of stars and 
gas in a lopsided potential, and reports that lopsided potentials arising from disks alone are not 
self-consistent; rather, a stationary lopsided disk must be responding to asymmetries in the halo. 
Finally, some numerical studies such as that done by Syer & Tremaine (1996) indicate that fluid 
disks may have non-axisymmetric equilibrium states, even when embedded in an axisymmetric 
halo. 

In certain galaxies, NGC 5474 for example (see van der Hulst & Huchtmeier 1979 and Rownd 
et al. 1994), the center of the optical light is separated from the center of the dynamical H I 
distribution, possibly indicating that galaxies do not rest in stationary states of the potential at 
all, but rather that the optical galaxy may be in a libration about the minimum of the the H I 
or dark matter potential. The preliminary numerical work reported by Levine and Sparke (1998) 
is suggestive that optical disks offset from dynamical centers can result in intriguing forms of 
non-axisymmetry. 

Several methods have been developed to quantify the frequency of asymmetry in galaxy 
images. Rix Sz Zaritsky (1995, hereafter RZ) and ZR analyzed the surface brightness distributions 
in if'-band images of spiral galaxies. Fourier signal strengths A m (r) for the m spiral modes were 
calculated, with large signal strengths (compared to the average surface brightness Aq) indicating 
asymmetry. RZ concentrated on the m = 1 "lopsided" Fourier mode as the primary indicator of 
asymmetry. Both RZ and ZR conclude that as many as 30% of galaxies show significant deviations 
from axisymmetry by this measure. 

Conselice (1997, hereafter C97) has proposed that asymmetries be measured by a method 
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in which galaxy images are rotated 180° and subtracted from the original image. Any residuals 
are interpreted as asymmetries. C97 does not provide error estimates of the derived asymmetry 
parameter, but comments that in the i?-band, at least 54% of 43 sample galaxies have asymmetry 
parameters A < 0.1, where A = and A = 1 are representative of maximal symmetry and 
asymmetry, respectively. The observation is also made that galaxies tend to be more asymmetric 
in the J-band than in the i?-band. 

Methods for measuring other than global asymmetries have also been proposed. For example, 
Elmegreen and Elmegreen (1995) study the symmetry of the two dominant inner spiral arms in a 
sample of 173 bright galaxies. In that study, the relative lengths and positions of the spiral arms 
were measured and compared. They find that in most galaxies, the inner spiral arm lengths are 
equal to within 20%, and are within 20° of being 180° apart. Such a method could be extended 
into the outer regions of a grand design spiral, in lieu of more global methods which can confuse 
strong spiral structure with asymmetry. 

In their study, BLS compared the H I surface density profiles on both sides of the galaxy for 
a sample of galaxies mapped in H I with sufficient resolution to measure the degree of asymmetry 
in H I disks. Using a much larger dataset of global H I profiles obtained from the literature, 
Richter & Sancisi (1994) classified asymmetry qualitatively. Based essentially on differences in the 
fluxes from the two horns, or total differences in flux between approaching and receding sides, 
they assigned galaxies to one of three categories of apparent asymmetry: "strong," "weak," or 
no asymmetry. Using the qualitative assignment, Richter & Sancisi concluded that only 47% 
± 5% of the global profiles show no asymmetry. Furthermore, they argued that an observed 
correspondence between obviously asymmetric H I maps (such as that for M101) and lopsided 
line profiles justifies the proposition that an observed asymmetry in the global H I profile can be 
taken as evidence for a non-circular density distribution. It should be noted, however, that this is 
a converse argument complicated by the fact that H I profiles contain both velocity field and H I 
distribution information, which are often difficult to disentangle. Richter h Sancisi did attempt 
to correct for the presence of companions or the presence of unusual kinematics by eliminating 
peculiar-looking profiles a priori. However, because the H I emission profiles used for this study 
were taken from observations made with a smorgasboard of telescopes, the effects of pointing error 
(which could manifest itself as an asymmetry in the profile) could not be determined. 

Haynes et al. (1998a) have conducted a similar study of H I line profiles in which spectra for a 
sample of 104 disk galaxies were taken using the 43-meter telescope at Green Bank and examined 
for asymmetry in a quantitative way. Using two methods for the determination of asymmetry, and 
after determining that pointing errors produced negligible deviations, the conclusion supported 
the previous results that asymmetry is found in ~ 50% of global H I profiles. However, those 
authors caution that asymmetry in line profiles obtained with single-dish telescopes can arise for 
a variety of reasons unrelated to disk asymmetry, such as contamination from companions within 
the main beam or sidelobes. 
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In this paper, we present a new photometric method for quantifying deviations from 
axisymmetry in optical images of disk galaxies. Essentially, we attempt to determine the likelihood 
that the light distribution in a galaxy is inherently symmetric by comparing the relative fluxes 
received in different directions from the optical center. In §2, we present the observational sample 
of 32 face-on and nearly face-on galaxies and describe the basic image acquisition and reduction 
procedures. In §3, the photometric method for the determination of departures from axisymmetry 
is introduced and then applied to the target sample. Finally in §4, the method is discussed in 
comparison with other quantitative determinations of asymmetry in spiral disks. 

2. Observations 

For the purpose of developing a method of measuring departures from axisymmetry, deep 
i?-band images were constructed for a set of 32 nearby face-on disk galaxies covering a range of 
spiral Hubble classes. In this section, we discuss the observational and image reduction procedures 
for the target sample of galaxies. 

2.1. Image Acquisition and Data Reduction 

The method of sectors described in the following sections was used to determine differences 
from symmetry for a sample of 32 face-on disk galaxies listed in Table || and displayed in Figure 
1. Observations were conducted with a Harris R-b&nd filter under nonphotometric conditions. 
Each galaxy except NGC 3393, NGC 3450, and MCG-5-34-002 was observed with the KPNO 0.9 
meter telescope configured to //7.5 and equipped with the T2KA direct imaging camera. The 
camera setup yielded a gain of 3.6 electrons per adu and a / .'68 per pixel scale. Only the inner 
1024 x 1024 pixels were illuminated, resulting in a field of view of 11'. 6 x 11'.6. Read noise in the 
camera was 4 electrons per pixel. Observations took place during several observing periods during 
November and December 1996 and February and March 1997. 

The three galaxies NGC 3393, NGC 3450, and MCG-5-34-002 were observed in April 1997 at 
the CTIOg 0.9 meter telescope with the T1K2 imaging camera. The gain in this camera is 3.35 
electrons per adu and the plate scale is / .'384 per pixel at Cassegrain focus. Read noise is 4.6 
electrons per pixel. 

At both telescopes, individual exposures of 600 seconds each were taken; the total integration 
time per object was varied depending on the atmospheric transparency. 
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Image reduction was accomplished in the /iMi^ image processing environment. These 
exposures were then debiased, sky-subtracted, and flatfielded using twilight flats. The exposures 
were scaled by measuring the integrated fluxes from 5-10 stars in the field. Because the 
transparency was variable, the frames were scaled to the single frame showing the greatest 
number of counts regardless of airmass. The frames were then combined using weights inversely 
proportional to the scale factor used, to produce images with total exposure times of between 
3000 and 4800 seconds. Because the images are not photometric (and in some cases, far from 
photometric), the relation between exposure times and total integration time is not clear. This 
variability does not affect the measurements conducted here, however, because the interesting 
quantities involve only magnitude differences of regions within individual images. Estimated errors 
which involve the integration time, however, may be affected. 

After combination, foreground stars were removed from the frame and replaced with a 
second-order polynomial interpolation of the surrounding background. A record of the number 
of pixels so replaced was kept, to be used as part of the final error estimate. Elliptical isophotes 
were then fit to the galaxy using a slightly modified version of the STSDAS^j ISOPHOTE package 
(Haynes et al. 1998b). The model galaxy thus obtained was used to determine each galaxy's center 
of light, inclination, and disk scale length R^. Scale length was determined by finding the best 
fit line in a surface brightness profile for the outer regions of a galaxy. Inclination was estimated 
under the assumption that galaxies seen face-on should appear circular. This assumption, 
although widely used to determine inclinations, is suspect in a project wherein inherent circular 
symmetry is not assumed. It is, however, justified statistically by studies such as that by Binney 
and de Vaucouleurs (1981) and in RZ, in which the observed distribution of galaxy ellipticities 
on the sky was used to estimate the distribution of the intrinsic axial ratios of disk and galaxies. 
We thus take the intrinsic ellipticities of disk galaxies as finite, but well within the errors of our 
ellipticity measurements. These values and the formal numerical errors are included in Table |l[ 
Following this basic data reduction, images were analyzed for asymmetry using the method of 
sectors described in the next section. 

The images shown in Figure [l] are individually scaled and trimmed for illustrative purposes, 
and most stars therein are masked to convey the visual appearance of each object. 

2.2. The Sample 

The observational sample selected for this study was chosen based on applicability for 
imaging for the purpose of employing the method of sectors discussed in The primary criterion 
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was face-on aspect, to eliminate confusion due to, for instance, the presence of dust lanes and 
extinction in the disk; the main source was the Third Reference Catalogue of Bright Galaxies 
(de Vaucouleurs et al. 1976; RC3). The galaxies in the observed sample were selected based 
on apparent disk size 1.'5 < D25 < 5', apparent ellipticity e* = R25 — 1 < 0.15, and narrow 
(Whi ^ 100 km s -1 ) HI line width. The latter restriction is a further diagnostic for low inclination 
in rotation-dominated disks. Since turbulence contributes to broadening at a level of ay (Broeils 
1992; Giovanelli et al. 1998) ~ 8 — 15 km s _1 and deviations from coplanar rotation may be quite 
common (Lewis 1987), H I line widths arising from spirals that are observed to be as narrow as 
100 km s _1 clearly suggest that the disk is seen face-on. Figure [2] demonstrates how rarely H I 
line widths are observed to be so narrow. The top panels shows the distribution of observed H I 
line widths W hi for a sample of 3123 galaxies in the Local Supercluster for which H I line widths 
are available either in the RC3 or, preferentially, in our own private database. Reported velocity 
widths have been converted to a systematic definition and corrected for instrumental broadening, 
signal-to-noise and redshift stretch, following Giovanelli et al. (1998) and Haynes et al. (1998c). 
The lower three panels show the distributions for spiral subtypes as indicated. Most of the narrow 
Whi profiles seen in the total sample arise from the low mass dwarf irregular galaxies, whereas 
only a small percent of spirals, those with inclinations < 15°, have narrow widths. We have thus 
used such widths as an important indicator of low inclination. 

Because of the desire to understand the nature of departures from symmetry in isolated disks, 
galaxies with no known nearby companions were preferentially included in the sample. In contrast, 
NGC 5474, which lies within the gravitational potential of M101 and is well-known for its high 
degree of asymmetry was included for comparative purposes. Likewise, NGC 1637 was added to 
the sample because of its well-known optical and infrared asymmetry despite its symmetric H I 
line profile (see Haynes et al. 1998a and Block et al. 1994). 

Because the observations were conducted during non-photometric conditions when the main 
program could not be carried out, final target selection was not controlled but rather dependent 
on the ill-fortune of the primary program. Table [l] lists relevant information about the program 
objects. Entries in the table til C c\S follows: 

Col. 1: Entry number in the Uppsala General Catalog (Nilson 1973), where applicable, or else in 
our private database, referred to as the Arecibo General Catalog (AGC). 

Col 2: NGC or IC designation, or other name, typically from the Catalog of Galaxies and Clusters 
of Galaxies (Zwicky et al. 1960), the ESO-Uppsala Catalog (Lauberts 1982) or the Morphological 
Catalog of Galaxies (Vorontsov-Velyaminov & Arhipova 1968). 

Cols. 3 and 4: Right Ascension and Declination in the 1950.0 epoch, typically from the AGC. In 
general, the listed positions have 15" accuracy. 

Col. 5: The morphological type code from the RC3. 

Col. 6: The major and minor diameters, D25 and 0I25, in arcmin, from the RC3. 
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Col. 7: The H I line width, Whi, corrected for instrumental broadening, smoothing, and 
signal-to-noise, in km s _1 . Values are taken from the literature as recorded in the AGC and have 
been converted with appropriate corrections to the system adopted by Haynes et al. (1998c) and 
Giovanelli et al. (1998). 

Col. 8: The inclination and its associated error, in degrees, derived from the ellipse fitting 
procedure described in §^l], under the assumption of an infinitely thin disk. The errors presented 
are formal numerical errors derived from constraints on geometrical quantities associated with 
the isophotes, including position angle and ellipticity. Because the galaxies are so nearly face-on, 
uncertainties in position angle are typically large. We therefore believe that the presented errors 
are overestimates. 

Col. 9: The heliocentric systemic velocity, V©, in km s" 1 , from the AGC. 

Col. 10: The total exposure time of the combined i?-band image. Note that the images were all 



some 



obtained under non-photometric conditions of differing transparency. As discussed in §2.1, 
compensation for the variation in conditions was made by varying the total exposure times and 
by weighting each exposure in forming the combination by the inverse of its scale factor. 

While the sample size is too small to treat its completeness, it was compared to the larger 
sample of all spiral galaxies of distance cz < 3000 km s ~ x in the RC3 with observed values of Z?25 
and i?25, with the subsample thereof including only those "face-on" galaxies with R25 < 1.3, and 
with the subsample comprised of those face-on galaxies within 2000 km s _1 . The distributions 
of several parameters were examined for each sample and mean values obtained. For a thorough 
investigation into the completeness of the catalogue, and distributions of parameters in the 
Local Supercluster, see Roberts & Haynes (1994). The parameters examined were: the apparent 
optical diameter £%, the apparent axial ratio R25, corrected bolometric magnitude B^, B — V&nd 
U — B colors, mean surface brightness m' e , heliocentric velocity, and type index. The analysis 
indicates that the sample galaxies, with (D25) = 3'.66 ± Of 9, are larger on the sky than on average 
for the supercluster ((-D25) = 2f 12 ± l'.O). The mean B — V&nd U—Bcolors and surface brightness 
of the sample are those of the local supercluster to within errors. The sample spans the spiral 
sequence with a mean type of Sc. While the sample is in no way complete, it is nonetheless 
representative of face-on galaxies within the Local Supercluster. 



3. The Method of Sectors 

We wish to construct a quantitative method for determining when a galaxy disk differs 
significantly from axisymmetry. We have decided on a method which is geometrically based, 
provides a global measure of asymmetry, and is applicable to a wide variety of data sets. 

The chosen method quantifies departures from disk axisymmetry by dividing the galaxy into 
a number of equal-area trapezoidal "wedges" in which photometry is to be performed (see Figure 
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H). The wedges are derived from a series of triangles, with the apex of each located at the center of 
light of the galaxy. The lateral boundaries of each wedge are radials emanating from this center. 
Each triangle is then truncated at a predetermined radius near the bulge, and extended outwards 
to the edge of the visible galaxy, thus creating a trapezoidal section. For most galaxies in the 
sample, inner and outer boundaries of Rd and 5Rd were chosen to ensure that the analysis would 
not be affected by the bulge component, especially including asymmetries introduced by bars, 
while including as much of the disk luminosity as possible. When the galaxy under consideration 
is slightly inclined, the sector pattern can be projected to the appropriate angle of inclination, 
generally 15° or 30°. Photometry is conducted in each of the sectors, and the magnitude of each 
sector is compared to that of the others. The magnitude difference between sectors i and j is 
AM%> = \Mi — Mj\, and the largest magnitude difference between sectors is AM™ 1 , where n 
represents the number of sectors used. AM™ 1 is the quantitative measure of asymmetry, and 
where this quantity is significantly different from zero, this is taken as a detection of such. 

To ensure that the method was sensitive to asymmetries in the visible flux, the method was 
tested by considering the very asymmetric galaxy NGC 5474 (Figure ||). This galaxy contains 
a nucleus which is so offset from the disk structure as to be disconnected from it. It does, 
however, have a roughly north-south line of approximate reflection symmetry. Using 6 sectors in 
a configuration aligned with the line of symmetry, we performed photometry in the sectors and 
determined the greatest difference between segments was AM™ 1 = 2™ 14 magnitudes. Because 
the sectors were aligned with the line of symmetry, we expected complimentary sectors to have 
similar total magnitudes. In this test, complimentary sectors across the line of symmetry differed 
only by an average of AMg =0^15 magnitudes. Similar results were obtained when the sector 
boundaries were allowed to vary away from perfect alignment with the symmetry axis by 5°, and 
for sets of 8 and 10 sectors, with average complimentary sector differences being AMg = CP20 
and AMf = CP23, respectively. The maximum magnitude differences in these cases were 
M max _ 2 m 34 and M m a x _ 2 m 36) consistent with AM 6 max to within the values of AM C . 

Errors in the measured flux of a region were estimated as the sum of two distinct factors. The 
first of these is the formal measurement error in the flux 5a including both the expected read noise 
and Gaussian sky noise. The second factor represents uncertainties in the masking of foreground 
stars. These two errors were added in quadrature for each galaxy to obtain an error estimate in 
the flux, and then transformed via the usual relations to a final error estimate in the magnitudes. 

The most ideal measure of the errors associated with this method would be to perform 
photometry on the sky in each frame, using the same sector pattern as used on the galaxy. The 
image frames, however, were generally not large enough to permit this. In order to understand 
the magnitude of the errors involved, therefore, the method of sectors was applied to a symmetric 
calibrator galaxy. For this calibration, it was important to choose a symmetric, flocculent disk 
galaxy so that neither spiral structure nor the presence of strong H II regions would contribute 
disproportionate amounts of luminosity to any one wedge. For these reasons, the galaxy 
MCG-1-60-011 was selected. Although classified as Sd in the RC3, this galaxy exhibits rather 
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tightly wound, very symmetric flocculent spiral structure and exhibits few H II regions. 

Observed values of the maximum flux difference Af max for MCG- 1-60-011 were interpreted as 
lex errors away from perfect symmetry in the fluxes. Furthermore, formal errors 5a were calculated 
for each flux determination. If /9 mas k is the density of masked pixels in the image, and A is the 
area of a sector, then the errors associated with masking should be proportional to Ap mas \ i . We 
therefore introduce the proportionality constant j3, such that when added in quadrature, the total 
error estimate becomes 



If the Af max values obtained for this galaxy are truly representative of symmetry, then we 
expect them to be of order a. Setting a = Af max a priori, it is possible to obtain a value of (3 for 
each segment size used. We find (3 to be roughly constant over segment size, thus the proposition 
of using Af max for the estimates of the errors is justified. A value of /3 = 92.3 ADU-pixel was 
taken for use with the remainder of the sample, as this is the average of the three values, none of 
which vary from the mean value by more than 5%. 

Table |2| contains the measured values of AM max for the 32 galaxies in the sample with la 
errors. Each galaxy was studied using three segment patterns consisting of 6, 8, and 10 sectors. 
Using a larger number of segments effectively compares smaller portions of the galaxy than when 
fewer segments are used, but in this regime, the estimated errors tend to increase. For the galaxies 
more than 5a away from symmetric, we have included the type of asymmetry seen based on the 
location of the minimum and maximum flux segments in the 8-segment runs. When the extreme 
segments are adjacent to one another, we desginate the galaxy "bisymmetric." When the extreme 
segments form a 135°or 180° angle, we call the galaxy "lopsided." Otherwise, the extreme segments 
form a 90° angle and the galaxy is "boxy." 

Of a total of 32 galaxies, 10 (31%) are asymmetric at the 5a level, while 19 (59%) are 
asymmetric at the 3a level, where "asymmetric" is taken to mean AM max > at the stated 
confidence level in at least one trial. Because the values of AM max are positive definite, however, 
it is important to note that any given measurement is more likely to be an overestimate than an 
underestimate of the "true" asymmetry. We therefore prefer to use the 5a values to report that 
approximately 30% of nearby field galaxies are optically asymmetric. This result agrees with the 
estimate found by RZ and ZR. 

Standard error estimates as described above could not be applied to NGC 5474, because the 
galaxy is so asymmetric that all attempts at isophotal ellipse fitting failed. Due to the similarities 
between the error estimates obtained with MCG- 1-60-011 and the AM C values obtained for 
NGC 5474, we adopted AM^ as the error estimates for that galaxy. If anything, these are probably 
overestimates. 




(1) 
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4. Discussion 

As illustrated in Table |2|, the degree of asymmetry returned by the method of sectors agrees 
roughly with a qualitative eyeball judgment of asymmetry for each galaxy. The only truly deviant 
result was obtained from NGC 1073, which while qualitatively only slightly asymmetric, returns 
a very significant (~ 28a) deviation from symmetry. We believe that this result is due the fact 
that although the galaxy appears relatively symmetric in the disk, the nucleus is offset from the 
center of the bar by about 10" to the northeast. It should be noted also that this galaxy exhibits 
a separate region of greatly enhanced surface brightness along the bar. Although this region did 
not formally affect the measurements, it emphasizes the fact that the bar region of NGC 1073 is 
peculiar. Such an offset of the nuclear region manifests itself as a strong lopsided asymmetry. 

In order to determine whether asymmetric galaxies share any other common characteristics, 
the subsample of the nine galaxies found to differ from symmetry by at least 5a (except for 
NGC 1073 and the peculiar galaxy NGC 5474) was compared to the overall sample of 32 galaxies. 
If we consider SAB galaxies to have .5 of a bar, then the average "barredness" of the asymmetric 
galaxies is about 0.44, compared to 0.53, the average "barredness" of the overall sample. We 
therefore find that lack of symmetry is about evenly distributed between barred and unbarred 
galaxies. The average morphological type of the asymmetric sample is Sbc. 

We also report on the dispositions of the five galaxies which were studied both by ZR 
and the current work. These galaxies are: NGC 600, NGC 991, NGC 1302, NGC 7742, and 
MCG-5-34-002 (also known as ESO 446 G 031). Of these five galaxies, RZ and ZR find that all 
but NGC 991 and MCG-5-34-002 have values of (Ai) < 0.2, which they interpret as indicative of 
overall symmetry. MCG-5-34-002 is reported to have a value (Ai) = 0.210 db 0.021. As stated in 
ZR, (A\) measurements are biased upwards, especially for small (Ai), so this measurement is not 
inconsistent with symmetry. These findings agree with the current work, in which only NGC 991 
is found to be different from symmetric by more than 5a. 

We also examine the galaxies in both our sample and the sample of C97. These galaxies are: 
NGC 3596, NGC 3631, NGC 4136, and NGC 5701, which in the latter study are given fi-band 
asymmetry parameters of 0.10, 0.12, 0.14, and 0.054 respectively. Although C97 does not provide 
a "cut-off" value for the asymmetry parameter A(R), we can see in that study that values range 
from A{R) = for the most symmetric galaxies to about A{R) = 0.16 for the most asymmetric. 
Our data are roughly consistent with that of C97, in that we find that the former three galaxies 
are asymmetric and the latter is symmetric. 

It is interesting to note also that global asymmetry may not necessarily be correlated with 
asymmetry in the inner regions, as measured by Elmegreen and Elmegreen (1995). Even the 
strongly asymmetric galaxy NGC 1637 is reported to have an inner arm length difference of only 
23%, and an arm position difference of only 20°, a result only slightly above the median in both 
cases. 
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The method of sectors has advantages over other methods in that it is neither dependent 
on assumed dominant even spiral modes, nor is it particularly sensitive to small deviations in 
inclination, which may manifest themselves in the magnitude of the m = 2 mode. This is not 
only because the sectors can be modified to account for slight inclinations, but also because even 
if the sectors are not modified, a slight deviation 5i in inclination will be detected as a slight 
decrease in the magnitude of sectors along the minor axis only of order 2.5 log (cos (Si)). This error 
is negligible if 2.5 log (sin (Si)) < 0™1 for most galaxies, requiring 5i <^ 25°. This method also, 
because of its purely geometric nature, is well suited to find asymmetries other than the m = 1 
"lopsided" type. Galaxies shaped like dumbbells, or with "boxy" or "triangular" shapes, which 
although asymmetric may be cut into equally luminous halves, would also be detected with this 
method. 

The method of sectors is suited both to the study of flocculent galaxies, and of galaxies with 
strong spiral arms, so long as those arms extend around the entire galaxy. This is because the 
integration over large sections tends to "wash out" effects of the spiral arms which might be seen 
as asymmetries in the Fourier method of RZ, for instance. We believe that this is a desirable 
effect, since we would not wish to detect symmetric spiral structure as an asymmetry. 

The method of sectors is also applicable to cases where the radially-dependant Fourier 
method is inapplicable or fails. The philosophy of the methods of RZ and C97 seems to be that 
the only truly symmetric galaxy is one where the surface brightness is a function of r alone. We 
believe that this definition ignores the central characteristic of disk galaxies, i.e. that they possess 
spiral arms, and even m = 1 spirals are not in general considered by observers to be qualitatively 
asymmetric. Thus these two methods would misclassify obviously symmetric, but predominantly 
one-armed spiral galaxies such as NGC 2326 (Figure ||) as lopsided. We therefore prefer to define 
symmetry based on equality of total integrated flux in various directions outward from the bulge. 
The method of sectors is designed to distinguish exactly such symmetry. 

Moreover, while the Fourier method is relatively easily applied to optical data from large 
nearby galaxies, it becomes progressively buried in noise as the data become more scarce. Such a 
method is not likely to be applicable to, for example: distant galaxies in the optical, HI synthesis 
data, or output of numerical modeling. The sector method does not suffer nearly as much from 
undersampling and allows direct comparison of optical photometry to other less prolific data sets. 
Often, numerical work is presented "as is," with few quantitative measures available to compare 
numerical results to actual data. In this regard, because of its use of flux integration over large 
regions, the method of sectors is easily adaptable to measurements of computer-generated data, 
allowing direct, quantitative comparison to natural phenomena. In the future, this method will be 
applied to aperture synthesis H I line mapping data and the results of numerical simulations of a 
subsample of the galaxies observed here. 

This research has been partially supported by NSF grant AST95-28860 to MPH, and 
AST93-20068 to RVEL, and has made use of the NASA/IPAC Extragalactic Database (NED) 
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a Disk inclination measured for this paper via isophotal ellipse fitting. (See text.) 
b Included in the sample of RZ. (See text.) 
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Table 2. Asymmetries Determined by Sector Comparison 
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indicates the detection of asymmetry in units of a. 

b "lopsided," "boxy," or "bisymmetric," if significantly asymmetric, (see text), 
indicates the qualitative appearance of asymmetry in a galaxy as "none," "slight," "moderate," 
or "strong." 

d Due to the extreme asymmetry of this galaxy, isophotal ellipses could not be fit. Errors 
expressed are values of AM", as described in the text. 
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Fig. 1. — Illustration of the 32 galaxies examined in the current paper and listed in Table 1. 
Relative sizes and intensities of the images are not to scale. 
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Fig. 2. — HI line width distributions for disk galaxies in the Local Supercluster, by morphological 
type. Note the small population density with HI line withs less than about 100 km s— 1. 

Fig. 3. — Illustration of the method used for quantifying asymmetry. Shown are galaxies 
NGC 1637 (left) and MCG-1-60-011 (right). The galaxy image is sectioned into some number 
of trapezoidal wedges centered on the nucleus in which photometry is performed. Some galaxies 
are not precicely face-on; in these cases, so long as the inclination is i ^ 30°, the sector pattern 
may be projected to the inclination of the galaxy, as is the case for NGC 1637 above. 

Fig. 4. — NGC 5474 with surrounding sector pattern. This galaxy, because of its extreme 
asymmetry, was used to test the sensitivity of the method of sectors. NGC 5474 lies within the 
gravitational field of M101, which lies about 1° to the NNW. 
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Fig. 5. — The relative surface brightness profile, in instrumental units, for the galaxies 
MCG-1-60-011 (top) and NGC 1637 (bottom). The dotted line in each panel shows the best 
linear fit to the outer portions of the profile. The derived disk scale length, Rd, is indicated by 
the dashed verticals. Each point is associated with an elliptical isophote. The average ellipticity of 
these isophotes over the region measured with the sector pattern is used to estimate the inclination 
of the galaxy. 



Fig. 6. — The galaxy NGC 2326 exhibits a single spiral arm in the outer regions which can be 
distinguished for three turns about the nucleus. Due to its strong m = 1 spiral mode, such a 
galaxy would be detected as lopsided with the methods of RZ and C97, but because its integrated 
luminocity is evenly distributed in 0, it is classified as symmetric by the present method. 
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